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Abstract 

The biological effect of one single radiation dose on a living tissue has been 
described by several radiobiological models. However, the fractionated ra- 
diotherapy requires to account for a new magnitude: time. In this paper we 
explore the biological consequences posed by the mathematical prolongation 
of a previous model to fractionated treatment. Nonextensive composition 
rules are introduced to obtain the survival fraction and equivalent physical 
dose in terms of a time dependent factor describing the tissue trend to- 
wards recovering its radioresistance (a kind of repair coefficient). Interesting 
(known and new) behaviors are described regarding the effectiveness of the 
treatment which is shown to be fundamentally bound to this factor. The 
continuous limit, applicable to brachytherapy, is also analyzed in the frame- 
work of nonextensive calculus. Also here a coefficient arises that rules the 
time behavior. All the results are discussed in terms of the clinical evidence 
and their major implications are highlighted. 
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1. Introduction 

As can be seen in jl| (and otlier works in the same issue) nonextensive 
Tsallis entropy ^ has become a successful tool to describe a vast class of 
natural systems. A recently developed model jsl of radiobiology shows this 
entropy definition could also be applied, not only to the development of living 
systems 0, [sf, but also to radiotherapy treatments. 

The new radiobiological model (maxent model in what follows) takes 
advantage of Tsallis entropy expression to describe the survival fraction as 
a functional of the radiation absorbed dose. This model is also based on a 
minimum number of statistical and biologically motivated hypotheses. 

The maxent model assumes the existence of a critical dose, Dq, that 
annihilates every single cell in the tissue. The radiation dose can be written 
as a dimensionless quantity in terms of that critical dose as x = (I/Dq, where 
d is the radiation dose. Then the support of the cell death probability density 
function, p{x), in terms of the absorbed dose x, becomes Q = [0; 1]. Tsallis 
entropy functional can be written. 



q- 



]f [x) dx 







(1) 



where q is the nonextensivity index. The survival fraction of cells will be 
given by / (x) = p (x) dx, that is the complement of the fraction of cells 
killed by radiation. In order to maximize functional ([1]) we must consider 
the normalization condition. 







p (x) dx = 1 (2) 

Also, following |6|, we must assume the existence of a finite g-mean value, 

-1 

(x) xdx = (3) 
Then the Lagrange multipliers method leads to, 

pix) = ^{i-xy-\ (4) 

with 7 = So, the survival fraction predicted by the model is 

f{x) = {l-x)\ (5) 
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valid for X G and requiring 7 > 1. 

This model has shown a remarkable agreement with experimental data 
js], 0], even in those limits where previous models are less accurate, mainly at 
high doses. The analysis of the model fit to experimental data also provides 
new hints about the tissue response to radiation: first, the interaction of a 
tissue with the radiation is universal and characterized by a single exponent 
(not dependent on the radiation type, energy or dose rate); second, the model 
includes a cutoff dose (this one, dependent on the characteristics of the ra- 
diation) above which every single cell dies. Furthermore, previous models 
can be derived as particular limiting cases. Finally, as for those models, its 
mathematical expression is simple and can be easily plotted and interpreted. 

The maxent model was derived for radiobiological survival fraction but its 
applicability could be extended to other processes. Indeed, every phenomena 
describable in terms of Tsallis entropy jsl, fulfilling the maximun entropy 
principle and exhibiting a critical cutoff (represented here by a; = 1), must 
follow (P. 

Nevertheless the expression 1^, understood as survival probability, lacks 
the extensivity property. In other words, for n events following the total 
survival probability should be found as a composition of the survival prob- 
abilities of the successive events. However, there is not a straightforward 
composition rule for those probabilities. 

Indeed, if two doses, xa and Xb are applied, the resulting probability from 
their composition has two possible values. If the total dose is assumed addi- 
tive, Jab = — xa — ■, that is, the individual probabilities under A and 
B events could not be treated as independent probabilities. Jab 7^ JaJb- On 
the other hand, if probabilities are multiplicative, / = (1 — xaY (1 — Xb)^ , 
doses would not fulfill the superposition principle for the equivalent physical 
dose, xab 7^ xa + Xb- 

The subject of this manuscript is to develop on the composition rules 
that would lead to the survival fraction and the equivalent physical dose of a 
fractionated processes, and to derive the biological implications of such rules. 
This will be approached within the frameworks of g-algebra and g-calculus 



9|, [lOl, lill, as far as they are the natural ones for the maxent model. 



2. Composition rules 

Each event described by ([5]) represents a measured energy impact, or dose, 
X causing an irreversible effect or hazard over a group of individual entities 
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leading to a survival probability of those entities. As (jSj) represents a nonex- 
tensive process, an appropiate set of composition rules must be developed in 
order to find the effect of several combined events on the group of entities. 

As it has just been exposed, if those composition rules are defined keeping 
the superposition principle for the dose, the probabilities are not indepen- 
dent of each other and vice versa, if the probabilities are multiplicative, the 
dose becomes non additive jl2|. Luckily, the nonextensive thermostatistics 



provides tools to find the right expressions in each case j9|, |lO|, lUl, Il2 |. 

If the survival probabilities are independent, the total probability for two 
events A and B is /ab = IaIb- So the nonextensive sum must be constructed 

a.sx®y = x + y — xy and we can write, 

XaB = Xa® Xb = Xa + Xb - XaXb ,ci\ 

_ (6) 

JAB — JAjB 

On the other hand, if the dose is additive, xab = xa + Xb, the nonextensive 
product must he x y = (x^/'^ + y^^'^ — l)"^ so we can write, 

Xab = xa + xb 

(7) 



The main issue here is that in clinical treatments both limits are not 
clearly distinct. Indeed, when events occur separate enough in time, tissue 
recovering capabilities make physical consequences of one of them indepen- 
dent from the others'. From a radiobiologist point of view this is similar to 
applying the next radiotherapy session after late effects of the former occur. 
However, if the events occur simultaneously the dose must be considered ad- 
ditive. In other words, ([6]) and ([7]) represent limit cases of the interaction 
process corresponding to t = oo and t = respectively, where t is the time 
between succesive events. 

In order to describe a real fractionated process, new generalized sum and 
product operators need to be introduced, taking into account that ([6]) and 
([7]) must hold in the multiplicative and additive limits, respectively. 

The resulting probability in ([6]) is the product of partial probabilities, and 
for the whole process, 

n 

F„ = n (1 - ^^v , (8) 

i=l 

where i runs along the events. 



4 



However, if the dose is additive, the total survival fraction follows. 



F„= (9) 



i=l 



Notice that it is possible to write ([2]) as a product, finding the expression 
that turns F after n — 1 events into F after n events. So, can be recast 
in the form, 

F-=(^- , Xi = n - , ^1 y (10) 

This expression can be interpreted as a modified ([S]) in which the denomina- 
tor, which plays the role of the annihilation cutoff, gets reduced, in practice, 
by an amount Xi after addition of the z-th event. On the other hand, for in- 
dependent events this critical cutoff would remain constant along the whole 
process. 

The new operators for nonextensive sum, and product, (^, must be 
defined to hold, 



subject to the condition ^^^i Xi — )■ Yli=i e — )■ 1- In this way, ([8]) and 

(fTOj) will be the limits of the new operators. Indeed, the coefficient e G [0, 1] 
acts as a session-coupling coefficient for equations (j8]) and (fTOjl such that 
e = 1 implies events are completely correlated while e = means they are 
fully independent, i.e. not coupled. 

Even though (fTTj) gives a closed and univocal definition of ^ and 
operators, this is an implicit definition. In order to use these operators an 
explicit definition is desired. 

The analytical expresion for the new operators © and ® can be found 
assuming there is a single event with an effective dimensionless dose X cor- 
responding to the whole process such that, 

F^ = {l-Xf = |l-0x, 1 (12) 

i=l 
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After the i-th event, the dimensionless effective dose would become, 



— ^i-l + Xi 



eX, 



i-l 



(13) 



assuming Xi = Xi. When the n-th event is given, then X„ = X. 
From this follows that, 



Xab = xa® xb = xa + xb 

fAB = Ia® Ib = Ia 



1 — X I 



1-tXA 

/ai/7_,(i_/^i/7) 



^ 7 



(14) 



and limit definitions (El) and ((71) are recovered for e = and e = 1 respectively. 
According to both limit interpretations, session-coupling e values will depend 
on the time between events and also on tissue repair or recovery capabilities. 



3. Biological and physical implications 

3.1. Isoeffect relationship 

One of the central concepts of radiotherapy is isoeffect relationships. An 
oncologist usually seeks treatments that produce the best outcome on the 
target tumor (Ftumor), while causing at most the maximum allowed damage 
on the surrounding healthy tissues {Fussue) ■ In other words, he seeks among 
the pairs of values (n, x) that give the same value of Fussue for the healthy 
tissue that one attaining the maximum value of Ftumor- Given the expression 
( IT2|) this can be reduced to find the pairs (n, x) that render the same value 
of X. 

Indeed, all fractionated treatments sharing the same value of effective 
dose, X, will provide the same value for the survival fraction. So, the same 
X will provide the isoeffect criterion for the fractionated therapy. 

In order to check the model reliability, it has been fitted to data from 



|13L |14| . |15| using a weighted least squares algorithm j3|. Those data sets are 
considered as a reliable source of clinical parameters (as the relation of 
LQ model |l6|). The results of the fit are shown in Figure [T] 

The obtained session-coupling coefficients show a survival fraction behav- 
ior far from the pure g-algebraic limits (e = and e = 1). Since e values for 
usual tissue reaction differ from limiting values, it is worth to further study 
the biophysical interpretation of this new parameter. 
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Figure 1: IsoefTect relationship data reported for mouse lung by (e = 0.50, -Do = 
11.3 Gy), mouse skin by [I^l (e = 0.58, Do = 24.0 Gy) and mouse jejunal crypt cells by 
[H (e = 0.62, Do = 16.1 Gy), fitted to (ini). 
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Every X value provides a different isoeffect relationsliip, as sliown in figure 
121 Once a treatment coefficient values (e and Dq) are known, the dosage can 
be tuned to obtain the desired effective dose by changing n and d. Notice 
that 7 does not play any role in this composition, thus reducing the number 
of model parameters to take into account here. 

As there is not enough experimental data available, in order to find 
session-coupling e values for known tissues or tumors we will use the LQ 
model of incomplete repairment to show how our model could be used to 
assess the desired therapy schedule. 

Let us suppose a healthy tissue H with 7 = 10.0 and Dq = 40.0 Gy 
surrounding a more resilient tumor T with 7 = 15.0 and Dq = 80.0 Gy. Now 
we will assume that H can not receive more than 36.0 Gy or X = 0.9. After 
finding the corresponding LQ model a and /3 values is easy to reproduce the 



isoeffect curves for incomplete repairment following [16| if the cell repair half 
time is known. We had chosen a repair half time of 3 hours for H and T 
but the same procedure could be applied for different repair half time values. 
Each of these curves represents a different treatment schedule characterized 
by the time (At) between sessions. From these curves the e values as a 
function of At could be found as shown in Figure [3l 

After the values of e have been determined for the tumor then the effective 
dose X received for each schedule could be found as shown in Figure HI This 
shows us that for small x values the best outcome is reached at more consec- 
utive sessions, whereas for more separated sessions the appropriate dosage 
is attained at higher x values. In particular, for the case of sample tissues 
H and T, described above, best results are found with a more fractionated 
treatment with its fractions scheduled as close as possible. 

Note that for a real example this procedure must be followed after finding 
the experimental values of e for each schedule. Even though illustrative, this 
example must be taken with caution as it is based on another model whose 
validity limits are not clear. 

3.2. Critical dosage 

Assuming the same physical dose per fraction, Xj = x, as is the case 
in many radiotherapy protocols, expression (fT3|) becomes that of a recur- 
sive map, describing the evolution of the effective dose in a treatment. The 
analysis of this map shows that, for every e there is a critical value of x, 

Xc = l -e, (15) 
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Figure 3: Session-coupling values found following the LQ model of incomplete repairment 
for an hypotetic tumor and healthy tissue as function of treatment time schedule. 
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Figure 4: X values as function of session adimensional dose x found for the hypotetic tumor 
T following the treatment schedules of incomplete repairment as function of treatment time 
schedule. Lines represents the approximate behaviour of X values. 
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Figure 5: The larger plot represents uq isolines as a function of x and e (dashed lines) 
above Xc{£) (solid line); below this line, killing all tissue cells is impossible. The small one 
represents critical values no in terms of Xc- 

dividing the plane (e, x) in two different regions (see figure[5]). For a treatment 
with X < Xc, there will always be a surviving portion of the tissue since 
always X„ < 1, for every n. However, if x > Xc, after enough fractions 
Xn > 1, meaning that effective dose has reached the critical value and every 
single cell of tissue has been removed by the treatment. Then it is possible 
to find no, the threshold value of n, that kills every cell, for a given therapy 
protocol. This is shown in the inset of Figure |5l 

If the desired result is the elimination of the radiated tissue cells, i.e. 
surrounding tissue is not a concern for treatment planning, Uq represents the 
minimum number of sessions needed to achieve this goal; any session after 
that will be unnecessary. On the contrary, if the therapy goal requires the 
conservation of tissue cells (for instance in order to preserve an organ), then 
the number of sessions must be lower than uq. 

The session-coupling parameter e is a cornerstone on isoeffect relation- 
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ships. A fractionated therapy of fully independent fractions requires a greater 
radiation dose per fraction, or more fractions, in order to reach the same iso- 
effect as a treatment with more correlated fractions. The session-coupling 
coefficient acts here as a relaxation term. Immediately after radiation dam- 
age occurs (e = 1) tissue begins to recover, as e decreases, until the tissue 
eventually reaches its initial radiation response capacity (e = 0). In other 
words, the formerly applied radiation results in a decrease of the annihilation 
dose (initially equal to Dq) describing the effect of the next fraction. The 
more coupled a session is to the previous one, the larger the value of e and, 
thus, the larger the effect on the critical dose will be. Notice that unlike 7, 
that characterizes the tissue primary response to radiation, e characterizes 
the tissue trend to recover its previous radioresistance. 

Correlation between fractions can be translated in terms of the late and 
acute tissue effects of radiobiology. Indeed, damaged tissue recovering capa- 
bilities should determine the value of e. Given a dosage protocol, to an early 
responding tissue would correspond e close to 0, whereas for a late respond- 
ing tissue, would be e closer to 1. Notice that in current working models for 
hyperfractionated therapies this repair and recovery effects are introduced as 
empirical correction factors [l7| . as will be required for the session-coupling 
coefficient. 

As it was shown in Q, nonextensivity properties of tissue response to 
radiation for single doses are more noticeable for higher doses than predicted 
by current models. On the contrary, for the same total dose, a lower dose per 
fraction will enhance nonextensive properties in fractionated therapies. In- 
deed, for high dosage a few fractions are applied in a treatment and a change 
in n is not required for different e values. However, in the lower dosage case, 
more radiation fractions need to be applied and the e parameter may be- 
come crucial. In this case n values move away from each other for isoeffect 
treatments with different e. So, in order to achieve the desired therapy ef- 
fects, fractionated radiotherapy must be planned for a tissue described by 7, 
varying x according to e. The session-coupling coefficient should be experi- 
mentally studied as its value tunes the annihilation dose along a radiotherapy 
protocol. 
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4. Continuous formulation 



4.1. Continuous limit 

For some radiation treatments as brachytherapy tlie irradiation is ap- 
plied in a single session but for a prolonged period of time. If the discrete 
irradiation sessions were close enough ( IT3l) could be written as, 



1 — X , , 

where r stands for the average absorbed radiation per unit time. At the early 
stages of continuous irradiation the effective dose is in general small, and is 
possible to assume eX < 1 and ~ 1 + eX. Then, 

X^r[l-(l-e)X], (17) 

where the terms of second order in eX and above have been neglected. How- 
ever, as can be seen in figure |6] this approximation moves away from (ITGll as 
time increases. 

4.2. Continuous irradiation 

It is obvious from dose additivity properties that in the continuous irra- 
diation case and for two time instants to and ti close enough, 

X= f\dt, (18) 

J to 

where r is the dose rate per unit time. However if both instants of time are 
far enough to make relevant the tissue recovering capabilities this expression 
becomes invalid. So, whereas a usual integration process could become valid 
in a short time period this is not true for longer intervals. So, in a similar way 
as was already done for the sum operation, a new definition for integration 
must be introduced. 

This can be done following jlo| and introducing the g-algebraic sum and 
difference, 

xSy = x + y — 9xy 
X B y 



(19) 

l^0y 



where 9 G [0, 1]. In those terms, a nonextensive derivative operation follows 
such that, 

dr t^to t-to I -Of ^ ' 
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Figure 6: Continuous limit approximation behavior for r = 0.1. As expected the solution 
of (fT6|) (thick continuous line) goes between a linear effect for e = 1 (thin continuous line) 
and the exponential approach to cutoff dose corresponding to e = (thin dashed line). 
Solution of ([T7)) is represented by the thick dashed line. 
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Then we can define the physical absorbed dose rate, r, as the nonextensive 
time derivative of the equivalent dose, 



Expression (12T|) can be rewritten as a standard ODE, 

X + OrX = r, (22) 

which can be solved in the usual way taking into account that 6 and r are 
in general functions of time. In the absence of recovering effects, the ap- 
plied effective dose would increase linearly, due to the applied radiation r. 
However a resistance force {6rX), that depends not only on tissue recovering 
characteristics but also on the dose rate and the effective dose itself, will slow 
down this increase. 

In order to illustrate the behaviour described by (l22l) . let us suppose r is 
constant (a common case in clinical practice) and 6 slowly varying in time, so 
that it can be also taken as a constant. Then it is straightforwardly obtained, 

X = l{l-exp(-0rt)}, (23) 

allowing to find the needed irradiation time to kill every cell in the tissue 
(^ = 1), 

and showing that effective dose increases at a decreasing speed, 

X = r exp {-9rt) , (25) 

until tissue cells get annihilated at time tk {X = 1). Under continuous 
irradiation, survival fraction decreases faster at the beginning of irradiation 
process. However, depending on dose rate and 6 coefficient, the killing process 
speed slows down until eventually every cell is killed. If the recovery capacity 
is very high [6 = 1] the radiation effects stack slowly and there will always 
be surviving tissue cells (tk = oo). Those radiation damages stack faster 
as long as tissue cells are less capable to recover themselves and if there 
is no recovery at all {6 = 0) the effective radiation dose grows linearly in 
time and cells get killed faster [tk = 1/t). This time shortening behavior 
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with decreasing recovering rate is also shown by other radiobiological models 



islljjj. 



Comparing fl22|) and f fT7|) we see that, in the limit of continuous dosage, 
they become the same expression with 6* ~ 1 — e. However this relation may 
become invalid at high exposures as effective dose becomes larger and eX 
becomes of order 1, as shown in figure El At this point, the fractionated and 
continuous treatments differ. 

This shows 6 could be considered constant only for a limited time of the 
continuous irradiation. It must be studied, in general, as a function of time, 
describing the growing resistance of tissue to be annihilated. This function 
should make that fl22p mimics the behavior of (fTB]) , shown in figure O 



5. Conclusions 

The use of Tsallis entropy and the maximum entropy ansatz (second law 
of thermodynamics) have allowed us to write a simple nonextensive expres- 
sion for the single dose survival fraction. The mathematical constraints, 
required to define the probabilities composition such that the two limiting 
behaviors are described, introduce a new parameter, relating the radiation 
sessions. The fits to available experimental data show that usual treatment 
have non trivial values of this parameter, i.e., are not close to the limit- 
ing behaviors. This makes the study of this coefficient relevant for clinical 
treatments and experimental setups. 

The existence of a varying critical dosage arises from these composition 
rules, providing a criterion to adjust the critical treatment that kills every 
tumor cell or minimize the damage caused to healthy tissue. This could 
be accomplished changing the number of sessions or the radiation dose by 
session, allowing to switch between isoeffective treatments. 

Also an expression for the effective dose in continuous irradiation treat- 
ments has been found, showing it is phenomenologically linked to the previous 
one. This has the potential to provide isoeffect relationships in continuous 
dose treatments such as brachytherapy. Besides, a relation between frac- 
tionated and continuous therapies could be established from the obtained 
coefficients. 
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